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Abstract 

The migration characteristics of carbon nanotubes from liquid phase to vapor phase in the refrigerant-based 
nanofluid pool boiling were investigated experimentally. Four types of carbon nanotubes with the outside 
diameters from 15 to 80 nm and the lengths from 1.5 to 10 pm were used in the experiments. The refrigerants 
include R1 13, R1 41 b and n-pentane. The oil concentration is from 0 to 10 wt.%, the heat flux is from 10 to 100 
kW-m" 2 , and the initial liquid-level height is from 1.3 to 3.4 cm. The experimental results indicate that the migration 
ratio of carbon nanotube increases with the increase of the outside diameter or the length of carbon nanotube. 
For the fixed type of carbon nanotube, the migration ratio decreases with the increase of the oil concentration or 
the heat flux, and increases with the increase of the initial liquid-level height. The migration ratio of carbon 
nanotube increases with the decrease of dynamic viscosity of refrigerant or the increase of liquid phase density of 
refrigerant. A model for predicting the migration ratio of carbon nanotubes in the refrigerant-based nanofluid pool 
boiling is proposed, and the predictions agree with 92% of the experimental data within a deviation of ±20%. 



Introduction 

Nowadays, the researchers show great interest in the pos- 
sible application of refrigerant-based nanofluids (i.e., the 
mixtures of nanopowders and conventional pure refriger- 
ants) for improving the performance of refrigeration 
systems. The researches showed that the refrigerant- 
based nanofluids have higher thermal conductivity than 
those of conventional pure refrigerants [1], the addition 
of nanoparticles enhances the solubility of mineral oil in 
HFC refrigerant [2], and the addition of nanoparticles 
can save energy consumption of air-conditioner and 
refrigerator [3,4] . Comparing with the spherical metal or 
metal oxide nanoparticles used in these researches, 
carbon nanotubes (CNTs) have one or two orders of 
magnitude higher in thermal conductivity, and CNTs can 
significantly enhance the thermal conductivity of base 
fluid [5-8] as well as the convective heat transfer coeffi- 
cient of base fluid [9], so CNTs have great potential for 
improving the performance of refrigeration systems. For 
applying CNTs in refrigeration systems, the phase- 
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change heat transfer characteristics of refrigerant-CNT 
nanofluid and the cycle behavior of CNTs in refrigeration 
systems should be known. The migration of CNTs from 
liquid phase to vapor phase in the pool boiling process of 
refrigerant-CNT nanofluid determines the distribution of 
CNTs concentration in the liquid phase and vapor phase, 
and then has significant effect on the phase-change heat 
transfer characteristics of refrigerant-CNT nanofluid as 
well as the cycle behavior of CNTs. Therefore, in order 
to evaluate the phase- change heat transfer characteristics 
of refrigerant-CNT nanofluid and the cycle behavior of 
CNTs, the migration of CNTs in the pool boiling process 
of refrigerant-CNT nanofluid should be researched. 

The migration of CNTs from liquid phase to vapor 
phase in the refrigerant-based nanofluid pool boiling 
can be divided to the following four physical processes: 

(1) the departure of bubble from the heating surface, 

(2) the movement of bubble and CNTs in the liquid 
phase, (3) the capture of CNTs by bubble, and (4) the 
escape of CNTs from the liquid-vapor interface. From 
the above analysis, it can be seen that the interaction 
between CNTs and bubble is the key factor causing the 
migration of CNTs from liquid phase to vapor phase. 
The existing flotation theory can accurately describe the 
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interaction between particles and bubbles. However, 
they can not directly used to predict the migration of 
CNTs from liquid phase to vapor phase during pool 
boiling because they are aimed at the conditions without 
phase change. The present investigation is beneficial to 
reveal the mass transfer mechanism of nano-scale solid 
powders from liquid phase to vapor phase during the 
phase-change process of fluid, and provides the theoreti- 
cal basis for evaluating the phase-change heat transfer 
characteristics of refrigerant-CNT nanofluid and the 
cycle behavior of CNTs. 

Until now, there is no published research on the 
migration characteristics of CNTs in the pool boiling of 
refrigerant-based nanofluid. The migration characteris- 
tics of nanopowders are mentioned only by the paper of 
Ding et al. [10], and are focused on one type of spheri- 
cal nanoparticle (CuO). In the paper, the authors found 
that the migrated mass of CuO nanoparticles in the 
pool boiling increases with the increase of the original 
mass of nanoparticles or the original mass of refrigerant. 
As the structure and the thermophysical properties of 
CNTs are different from those of spherical nanoparti- 
cles, the migration characteristics of CNTs in the pool 
boiling of refrigerant-based nanofluid may be different 
from those of spherical nanoparticles, and should be 
investigated. 

The existing researches on the pool boiling heat trans- 
fer characteristics of nanofluids containing CNTs can be 
divided into three categories as follows: (1) Pool boiling 
heat transfer of water-CNTs nanofluids. The experimen- 
tal results showed that CNTs can enhance the pool boil- 
ing heat transfer of water [11-13] or deteriorate the pool 
boiling heat transfer of water [14], and the influence of 
CNTs on the pool boiling is related to CNTs concentra- 
tion [13]. (2) Pool boiling heat transfer of pure refriger- 
ant-CNTs nanofluids. The experiments by Park and 
Jung [11,15] showed that CNTs can enhance the pool 
boiling heat transfer of pure refrigerants (R22, R123, 
and R134a), and the enhancement is related to the heat 
flux. (3) Pool boiling heat transfer of refrigerant-oil mix- 
tures with CNTs. Experiments by Peng et al. [16] 
showed that CNTs can enhance the pool boiling heat 
transfer of refrigerant-oil mixtures, and the enhance- 
ment increases with the decrease of CNTs' outside dia- 
meter or CNTs' nanolubricant mass fraction, while 
increases with the increase of CNTs' length or CNTs' 
mass fraction in the CNTs' nanolubricant. From the 
above researches, it can be seen that the CNTs physical 
dimension (i.e., the outside diameter and the length of 
CNTs), refrigerant type, CNTs concentration, oil con- 
centration, and heat flux have influences on the pool 
boiling heat transfer of refrigerant-based nanofluid. 
Therefore, the influences of the above factors on the 
migration characteristics of CNTs need be concerned. 



In addition, the initial liquid-level height affects the pool 
boiling heat transfer, so the influence of initial liquid- 
level height on the migration characteristics of CNTs 
also needs be concerned. 

The objective of this paper is to experimentally inves- 
tigate the influences of CNTs physical dimension, refrig- 
erant type, oil concentration, heat flux, and initial 
liquid-level height on the migration characteristics of 
CNTs in the refrigerant-based nanofluid pool boiling at 
different original CNTs concentrations, and to propose 
a model for predicting the migration ratio of CNTs in 
the refrigerant-based nanofluid pool boiling. 

Experiments 

Test conditions and experimental objects 

Test conditions are divided into five categories, as tabu- 
lated in Table 1. 

The objective of category 1 is to investigate the influ- 
ences of CNTs physical dimension on the migration 
characteristics of CNTs. Four types of CNTs with differ- 
ent physical dimensions (numbered as CNT#1, CNT#2, 
CNT#3, and CNT#4) produced by the chemical vapor 
deposition method are used in these test conditions. 
The physical dimensions of these four types of CNTs 
are shown in Table 2 and the TEM (transmission elec- 
tron microscope) photographs of the CNTs are shown 
in Figure 1. In these test conditions, the other influence 
factors including the refrigerant type, oil concentration, 
heat flux and initial liquid-level height are fixed. 

The objective of category 2 is to investigate the influ- 
ences of refrigerant type on the migration characteristics 
of CNTs. Three types of refrigerants including R113, 
R141b, and n-pentane are used in these test conditions, 
belonging to CFC refrigerant, HCFC refrigerant, and 
alkane refrigerant, respectively. The reasons for choosing 
these three types of refrigerants are as follows: (1) R113, 
R141b and n-pentane are in liquid state at room tem- 
perature and atmospheric pressure while the widely 
used refrigerants (e.g., R410A) are in vapor state, so it is 
much easier to prepare refrigerant-based nanofluids 
based on R113, R141b, or n-pentane. (2) These three 
types of refrigerants have different chemical and ther- 
mophysical properties including molecular mass, density, 
dynamic viscosity, etc. The properties of these three 
refrigerants are given in Table 3. In these test condi- 
tions, the other influence factors including the CNTs 
physical dimension, oil concentration, heat flux and 
initial liquid-level height are fixed. 

The objective of category 3 is to investigate the influ- 
ences of oil concentration on the migration characteris- 
tics of CNTs. The lubricating oil RB68EP is used in the 
experiments. RB68EP is an ester oil with a density of 
0.964 g-cm" 3 at 15°C and kinematic viscosities of 66.79 
and 8.23 mm 2 's _1 at 40°C and 100°C, respectively, as 
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Table 1 Test conditions 



Objective of investigation 


CNTs 
type 


Refrigerant 
type 


Oil concentration 

x Q (wt.%) 


Heat flux q 
(kWrn 2 ) 


Initial liquid-level 
height L (cm) 


Original CNTs 
concentration 4> n (vol.%) 


Influence of CNTs physical 
dimension on migration 


CNT#1, 
CNT#2, 
CNT#3, 
CNT#4 


Rl 13 


0 


50 


2.0 


0.56, 1.11, 1.65, 2.19, 2.72, 
3.25, 3.77 


Influence of refrigerant type on 
migration 


CNT#2 


R113, 

R141b, 

n-pentane 


0 


50 


2.0 


0.56, 1.11, 1.65, 2.19, 2.72, 
3.25, 3.77 


Influence of oil concentration on 
migration 


CNT#2 


Rl 13 


1, 5, 10 


50 


2.0 


0.56, 1.11, 1.65, 2.19, 2.72, 
3.25, 3.77 


Influence of heat flux on 
migration 


CNT#2 


Rl 1 3 


0 


10, 20, 
50, 100 


2.0 


0.56, 1.11, 1.65, 2.19, 2.72, 
3.25, 3.77 


Influence of initial liquid-level 
height on migration 


CNT#2 


Rl 13 


0 


50 


1.3, 2.0 
2.7, 3.4 


0.56, 1.11, 1.65, 2.19, 2.72, 
3.25, 3.77 



reported by the manufacturer. The oil concentration is 
from 0 to 10 wt.%, covering the oil concentration in the 
actual refrigeration system. In these test conditions, the 
other influence factors including the CNTs physical 
dimension, refrigerant type, heat flux and initial liquid- 
level height are fixed. 

The objective of category 4 is to investigate the influ- 
ences of heat flux on the migration characteristics of 
CNTs. The heat flux is from 10 to 100 kW-m" 2 , covering 
the heat flux in the actual refrigeration system. In these 
test conditions, the other influence factors including the 
CNTs physical dimension, refrigerant type, oil concen- 
tration and initial liquid-level height are fixed. 

The objective of category 5 is to investigate the influ- 
ences of initial liquid-level height on the migration char- 
acteristics of CNTs. The initial liquid-level height is from 
1.3 to 3.4 cm. In these test conditions, the other influence 
factors including the CNTs physical dimension, refriger- 
ant type, oil concentration, and heat flux are fixed. 

Experimental apparatus 

The experimental apparatus used for testing the migra- 
tion characteristics of CNTs in the refrigerant-based 
nanofluid pool boiling mainly consists of a pool boiling 
device, a capture cover and a digital electronic balance, 
as schematically shown in Figure 2. The pool boiling 
device mainly consists of a boiling vessel and an electric 
heating membrane. The boiling vessel is a cylindrical 
glass container with the inside diameter of 50 mm and 
the height of 95 mm. The vessel is insulated with glass 



fibers to reduce heat loss to the surroundings. The elec- 
tric heating membrane is connected with the direct- 
current voltage power supply. The ampere meter with 
the calibrated precision of 0.5% is used for reading elec- 
tric current supplied to the heating surface, and a data 
acquisition system with the calibrated precision of 
0.002% is used to measure the voltage across the heat- 
ing surface. The heat flux through the heating surface is 
controlled by adjusting the heating power of the plate 
heater, and is calculated by the measured electric cur- 
rent, voltage, and heating surface area. The uncertainty 
of heat flux is estimated to be smaller than 1.2%. The 
capture cover is used to collect the CNTs spouted to 
the environment. The measurement range of the digital 
electronic balance is from 10.0 mg to 210.0000 g, and 
the maximum error is 0.1 mg. All the experiments are 
performed at atmospheric pressure (101.3 kPa) by vent- 
ing the boiling vessel to ambient. 

Experimental method 

The objective of the measurements is to get the 
migrated mass of CNTs from liquid phase to vapor 
phase in the refrigerant-based nanofluid pool boiling. 
Ding et al. [10] have proposed the weighing method to 
obtain the migrated mass of spherical nanoparticles, and 
this method is also used in the present study to get the 
migrated mass of CNTs. 

The experimental procedure for the refrigerant-based 
nanofluid without oil consists of the following steps: 
(1) adding CNTs with the original mass of m n0 to the 



Table 2 Physical dimensions of CNTs 

Property Mean outside diameter (cf out ) Mean inside diameter [d in ) Mean length Aspect ratio 

(/) (//Cfout ) 

CNT#1 15 nm 10 nm 1.5 Mm 100 

CNT#2 15 nm 10 nm 10 pm 666.7 

CNT#3 80 nm 18 nm 1.5 Mm 18.8 

CNT#4 80 nm 18 nm 10 pm 125 
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200 nm 




Si 

(a) CNT#1, J out =15nm, /=1.5|xm 




















200 nm 



(b) CNT#2, d out =15nm, /=10|xm 
200 nm 




200 nm 



(c) CNT#3, d out =80nm, 1=1 .5pm (d) CNT#4, d out =80nm, /=10|im 

Figure 1 TEM photographs of CNTs. (a) CNT#1; (b) CNT#2; (c) CNT#3; (d) CNT#4. 



boiling vessel; (2) weighing the total mass of the boiling 
vessel and the CNTs, (3) adding refrigerant with the 
original mass of m r0 to the boiling vessel; (4) opening 
the direct-current voltage power supply and heating the 
refrigerant-based nanofluid to be boiling; (5) adjusting 
the voltage to control the heat flux; (6) weighing the 
total mass of the boiling vessel and the CNTs, m 2 , when 
the refrigerant is entirely evaporated (the signal for the 
entire evaporation is that the mass of mixture does not 
change for 12 h); (7) calculating the migrated mass of 
CNTs with the equation Am n = mi - m 2 . 

The experimental procedure for the refrigerant-based 
nanofluid with oil consists of the following steps: (1) add- 
ing CNTs with the original mass of m n0 and oil with the 
mass of m Q to the boiling vessel; (2) weighing the total 
mass of the boiling vessel, the CNTs and oil, m 3 ; 



(3) adding refrigerant with the original mass of m r0 to 
the boiling vessel; (4) opening the direct-current voltage 
power supply and heating the refrigerant-based nano- 
fluid with lubricating oil to be boiling; (5) adjusting the 
voltage to control the heat flux; (6) weighing the total 
mass of the boiling vessel, the CNTs and oil, m 4 , when 
the refrigerant is entirely evaporated (the signal for the 
entire evaporation is that the mass of mixture does not 
change for 12 h); (7) calculating the migrated mass of 
CNTs with the equation Am n = m 3 - m 4 . 

Data reduction and uncertainty 

In order to quantitatively evaluate the migration degree of 
CNTs, the migration ratio of CNTs, & is defined by Eq. 1: 



? = Am n /m n0 



(1) 



Table 3 Properties of refrigerants in the experiments 



Property 


R113 


R141b 


n-Pentane 


Chemical formula 


CI2FC-CCIF2 


CH3CCI2F 


C5H12 


Molecular mass (g/mol) 


187.37 


116.95 


72.15 


Normal boiling point (°C) 


47.6 


32.05 


36.1 


Liquid-phase density (kg-m -3 ) 


1508 


1220 


606 


Vapor-phase density(kg-m" 3 ) 


7.4 


4.8 


3.2 


Liquid-phase dynamic viscosity (Pa-s) 


5 x 10" 4 


3.78 x 10" 4 


1.91 x 10" 4 


Liquid-phase isobaric specific heat (J-kg" 1 -K _1 ) 


940 


1,165 


2,441 


Surface tension (N-m -1 ) 


0.0147 


0.0173 


0.0138 
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Figure 2 Schematic diagram of experimental apparatus 



where, Am n is the migrated mass of CNTs, and m n0 is 
the original mass of CNTs. 

The original CNTs concentration is defined as the ori- 
ginal volume fraction of CNTs in the liquid phase 
(liquid refrigerant or liquid refrigerant-oil mixture), pre- 
sented as Eq. 2: 



Wno/Pn 



m n0 /p n + m r0 /p r ,L + m 0 /p 0 



(2) 



where, m r0 , m Q are the original mass of refrigerant and 
the mass of oil, respectively; p n , p r> L , p G are the density 
of CNTs, liquid-phase refrigerant and oil, respectively. 

The oil concentration is defined as the original mass 
fraction of oil in the liquid refrigerant-oil mixture, pre- 
sented as Eq. 3: 



m r0 + ra 0 



(3) 



The relative uncertainty of migration ratio of CNTs is 
calculated as: 



8 1 



——) 8Am n 2 + ( — 
Am n / \ m no 



<5ra n0 2 



(4) 



Determined by the accuracy of the digital electronic 
balance, the maximum uncertainties of the measured 
migrated mass of CNTs (SAm n ) and original mass of 



CNTs (Sm 0 ) are 0.2 mg and 0.1 mg, respectively. The 
maximum relative uncertainty of migration ratio of 
CNTs is obtained at the condition of the smallest 
migrated mass of CNTs and the migrated mass of 
CNTs, and calculated to be 2.5%. 

Tests under several conditions were repeated for three 
times, and it shows that the differences among the three 
testing results under each condition are less than 3%. 
Therefore, the experimental results are reproducible. 

Experimental results and analysis 

Influence of CNTs physical dimension on the migration of 
CNTs 

Figure 3 shows the migration ratio (£) of CNTs as a 
function of original CNTs concentration (0 n ) for CNTs 
with different physical dimensions. The values of C 
under these test conditions are in the range of 5.1% to 
approximately 27.8%. For fixed CNTs physical dimen- 
sion, C decreases with the increase of <fi n . For example, C 
for CNT#1 decreases by 15.8% with the increase of <p n 
from 0.56 to 3.77 vol.%. 

From Figure 3 it can be seen that the migration ratio 
(0 of CNTs increases with the increase of the outside 
diameter of CNTs (d out ) when the length of CNTs (/) is 
fixed. For example, at the condition of / = 1.5 urn, the 
value of C increases by maximally 286.4% with the 
increase of d out from 15 to 80 nm. It also can be seen 
from Figure 3 that the migration ratio (Q of CNTs 
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Figure 3 Influence of CNTs physical dimension on the 
migration ratio of CNTs. 
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Figure 4 Influence of refrigerant type on the migration ratio of 
CNTs. 



increases with the increase of the length of CNTs (/) 
when the outside diameter of CNTs (d out ) is fixed. For 
example, at the condition of d out = 15 nm, the value of 
C increases by maximally 25.7% with the increase of / 
from 1.5 to 10 [im. The possible reasons for the above 
phenomenon are as follows: The capture of CNTs by 
the bubbles generated in the pool boiling leads to the 
migration of CNTs. Brownian diffusion, interception, 
gravity settling and inertial impaction are four mechan- 
isms for the capture of particles by bubbles [17]. As the 
CNTs do not exhibit Brownian motion due to their high 
aspect ratio [18,19], the capture efficiency of CNTs by 
bubbles can be considered as the sum of the capture 
efficiencies caused by interception, gravity settling, and 
inertial impaction. Each of the above three captured effi- 
ciencies increases with Stokes diameter of CNTs. The 
increase of the outside diameter or the length of CNTs 
causes the increase of Stokes diameter of CNTs [20], 
thus the captured efficiency of CNTs by bubbles 
increases, which leads to the migration ratio (0 of 
CNTs increasing with increase of the outside diameter 
or the length of CNTs. 

Influence of refrigerant type on the migration of CNTs 

Figure 4 shows the migration ratio (0 of CNTs as a 
function of original CNTs concentration (0 n ) for differ- 
ent types of refrigerants. The values of C under these 
test conditions are in the range of 3.8% to approximately 
8.2%. For fixed refrigerant type, f decreases with the 
increase of (j) n . For example, C for R141b decreases by 
19.5% with the increase of <p n from 0.56 to 3.77 vol.%. 



From Figure 4 it can be seen that the migration ratio 
(0 of CNTs are in the order of R141b > R113 >n-pen- 
tane, the value of C in the R141b-based nanofluid is by 
maximally 10.7% larger than that in the R113-based 
nanofluid, and is by maximally 77.4% larger than that in 
n-pentane-based nanofluid. The possible reasons for the 
above phenomenon are as follows: (1) The dynamic 
viscosity values for these three refrigerants are in the 
order of R113 (5 x 10" 4 Pa-s) > R141b (3.78 x 10~ 4 Pa-s) 
>n-pentane (1.91 x 10" 4 Pa-s). The larger dynamic visc- 
osity causes the smaller capture efficiencies caused by 
gravity settling and inertial impaction, which leads to 
the smaller migration ratio of CNTs. (2) The liquid- 
phase density values of for these three refrigerants are 
in the order of R113 (1,508 kg-m" 3 ) > R141b (1,220 
kg-m" 3 ) >n-pentane (606 kg-m" 3 ). The larger liquid-phase 
density means the larger mass of liquid-phase refrigerant 
at fixed liquid-level height, thus the amount of bubbles 
generated in the pool boiling is larger, which leads to 
the larger migration ratio of CNTs. The influence of 
refrigerant type on the migration ratio of CNTs is deter- 
mined by the conjunct role of the above two aspects, 
and follows the order of R141b > R113 >n-pentane. It 
can be concluded that the migration ratio of carbon 
nanotube increases with the decrease of dynamic viscos- 
ity of refrigerant or the increase of liquid-phase density 
of refrigerant. 

Influence of oil concentration on the migration of CNTs 

Figure 5 shows the migration ratio (0 of CNTs as a 
function of original CNTs concentration (cj) n ) for 
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o I " " 1 J " 1 " " 1 

0.0 0.5 LO 1.5 2.0 2.5 3.0 3.5 4.0 4.5 

P n (%) 

Figure 5 Influence of oil concentration on the migration ratio 
of CNTs. 



different oil concentrations (x Q ). The values of £ under 
these test conditions are in the range of 1.3% to 
approximately 7.4%. For fixed x of £ decreases with the 
increase of (j) n . For example, at the condition of x 0 = 5 
wt.%, £ decreases by 28.2% with the increase of (j) n from 
0.56 to 3.77 vol.%. 

From Figure 5 it can be seen that the migration ratio 
{£) of CNTs decreases by maximally 70.7% with the 
increase of oil concentration (x Q ) from 1 to 10 wt.%. 
The possible reasons are as follows: (1) The dynamic 
viscosity and surface tension of lubricating oil RB68EP 
are larger than those of pure refrigerant, causing the 
dynamic viscosity and surface tension of liquid-phase 
refrigerant-oil mixture increase with the increase of x Q . 

(2) The increase of dynamic viscosity of liquid-phase 
refrigerant-oil mixture results in the decrease of capture 
efficiencies caused by gravity settling and inertial impac- 
tion, which leads to £ decreasing with the increase of x Q . 

(3) The increase of surface tension of liquid-phase 
refrigerant-oil mixture causes the increase of bubble 
departure diameter in the pool boiling, thus the capture 
efficiencies caused by interception, gravity settling, and 
inertial impaction decrease, which leads to £ decreasing 
with the increase of x Q . 

Influence of heat flux on the migration of CNTs 

Figure 6 shows the migration ratio (£) of CNTs as a 
function of original CNTs concentration (0 n ) for differ- 
ent heat fluxes {q). The values of £ under these test con- 
ditions are in the range of 5.5% to approximately 9.2%. 
For fixed q, £ decreases with the increase of 0 n . For 




4 I i i i i i i i i I 

0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 

% {%) 

Figure 6 Influence of heat flux on the migration ratio of CNTs. 



example, at the condition of q = 10 kWm" 2 , £ decreases 
by 21.5% with the increase of (j) n from 0.56 to 3.77 vol.%. 

From Figure 6 it can be seen that the migration ratio 
{£) of CNTs decreases by maximally 33.9% with the 
increase of heat flux (q) from 10 to 100 kWm" 2 . The 
possible reasons are as follows: (1) The increase of heat 
flux causes the increase of the velocity of departure bub- 
ble [21], thus the velocity of rising bubble increases. 
(2) The increase of velocity of rising bubble results in 
the decrease of capture efficiency caused by gravity set- 
tling, which leads to the decrease of £ (3) The increase 
of velocity of rising bubble results in the increase of 
capture efficiency caused by inertial impaction, which 
leads to the increase of £ (4) The increase of velocity of 
rising bubble results in the decrease of the bubble rising 
time in the liquid phase, causing the decrease of the 
amount of CNTs captured by bubbles, which leads to 
the decrease of £ The conjunct role of the above aspects 
leads to the migration ratio of CNTs {£) decreasing with 
the increase of heat flux. 

Influence of liquid-level height on the migration of CNTs 

Figure 7 shows the migration ratio (£) of CNTs as a 
function of original CNTs concentration (0 n ) for differ- 
ent initial liquid-level heights (L). The values of £ under 
these test conditions are in the range of 3.2% to 
approximately 19.8%. For fixed L, £ decreases with the 
increase of 0 n . For example, at the condition of L = 
2.7 cm, £ decreases by 16.5% with the increase of 0 n 
from 0.56 to 3.77 vol.%. 

From Figure 7 it can be seen that the migration ratio 
(£) of CNTs increases by maximally 446.9% with the 
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Figure 7 Influence of initial liquid-level height on the 
migration ratio of CNTs. 



increase of initial liquid-level height (L) from 1.3 to 3.4 
cm. The possible reasons are as follows: (1) The increase 
of initial liquid-level height causes the increase of the 
bubble rising time in the liquid phase, thus the amount 
of CNTs captured by bubbles, which leads to C increas- 
ing with the increase of liquid-level height. (2) The 
increase of initial liquid-level height causes the increase 
of CNTs escape probability on the liquid-vapor inter- 
face, which leads to C increasing with the increase of 
liquid-level height. 

Prediction of migration ratio of CNTs in the 
refrigerant-based nanofluid pool boiling 

As there is no published literature on the model for pre- 
dicting the migration ratio of CNTs in the refrigerant- 
based nanofluid pool boiling, the development of a new 
model is needed. The CNTs physical dimension, refrig- 
erant type, oil concentration, heat flux, and initial 
liquid-level height are five important factors influencing 
the migration of CNTs, and should be reflected in the 
new model. 

From the beginning of the pool boiling t 0 to the 
moment of t, the migration ratio of CNTs can be 
expressed as: 



m 



n,t 

m n0 



(5) 



where, m n>t is the mass of CNTs in the liquid phase at 
the moment of t; m n0 is the original mass of CNTs. 

According to the principle of mass conservation, the 
mass of CNTs in the liquid phase changed with time 



can be expressed as: 
dm n , t 



dt 



-Kn b m nit 



(6) 



where, K is the migration coefficient of CNTs (i.e., the 
migration proportion of CNTs caused by single bubble); 
n b is the number of generated bubbles per unit time. 

From Eq. 6, the following equation can be obtained. 



m n0 • exp (-Kn b t) 



(7) 



Combing Eqs. 5 and 7, ft 0 -t can be calculated as: 

f to _ £ = 1 - exp (-Kn b t) (8) 

Therefore, from the beginning to the end of pool boil- 
ing, the migration ratio of CNTs is: 



? = 1 - exp (-KN b ) 



(9) 



In Eq. 9, N b is the total number of generated bubbles 
from the beginning to the end of pool boiling. 

1. The calculation of K 

The migration of CNTs can be considered as the cap- 
ture of CNTs by bubbles combining the escape of CNTs 
from the liquid-vapor interface. In order to describe the 
capture process of CNTs by bubbles, the capture effi- 
ciencies of CNTs caused by interception, gravity settling 
and inertial impaction should be included in K. In order 
to describe the escape process of CNTs from the liquid- 
vapor interface, the escape probability of CNTs should 
be included in K. The original CNTs concentration has 
influences on the bubble diameter and bubble rising 
velocity during the pool boiling process of refrigerant- 
based nanofluid, and then has influence on the migra- 
tion of CNTs. Therefore, the CNTs concentration 
impact factor should also be included in K. The expres- 
sion of K is as follows: 



K = (oti + a G + oti N )Py 



(10) 



where, a h a G , and ai N are the capture efficiencies of 
CNTs caused by interception, gravity settling, and iner- 
tial impaction, respectively; /? is the escape probability of 
CNTs; y is the CNTs concentration impact factor. a h 
a G , a IN , /3, and / are expressed as Eqs. 11 to 15, respec- 
tively. 



a\ = a,2 



g(Pn~ PL) d S 



(11) 



(12) 
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q?in = #3 • St 



D 



1 #3 



ds 2 p n Ub 
18dbML 



H-L/2 



y = a 4 • <Pn 



(13) 



(14) 



(15) 



In Eqs. 11 to 15, D and are the bottom diameter 
and height of boiling vessel, respectively; p L , p b , p n are 
the density of liquid refrigerant-oil mixture, bubble and 
CNT, respectively; ft L is the dynamic viscosity of liquid 
refrigerant-oil mixture; d s is the Stokes diameter; d b is 
the diameter of bubble; u h is the bubble rising velocity. 
The calculation of p L and m L are presented in -Table 4. 

d s can be calculated by Eq. 16 [20]: 



U S — u out 



(16) 



d h can be calculated by Cole-Rohsenow correlation 
[22], as shown in Eq. 17: 



d b = 4.65 x 10 



-4j-PLCp,lTsat^ A 



5 

")4 



- Pb) 



(17) 



u h can be calculated by Malenkov correlation [21], as 
shown in Eq. 18: 



lag 



I 1 / 2 



4g(PL ~ Pb) ^ 

2(PL + Pb) 4(PL + Pb)J frfgPb 



(18) 



In Eqs. 17 and 18, C p> L , is the isobaric specific heat of 
liquid refrigerant- oil mixture; a is the surface tension of 
refrigerant-oil mixture; T SSLt is the saturation tempera- 
ture; hf g is the latent heat of vaporization. The calcula- 
tion of C Pt l and a are presented in -Table 4. 

2. The calculation of N b 

As the original mass of refrigerant is equal to the total 
mass of generated bubbles from the beginning to the 
end of pool boiling, N h can be calculated by Eq. 19. 



N b 



3D 2 H{l-x 0 )p L 
2pt,db 3 



(19) 



The nine coefficients of a±t ^3, ^4, bi, b^ b%, b^ 
and b 5 in Eqs. 11 to 15 are fitted based on total 105 
experimental data in this study. By nonlinear program- 
ming; solution method, the nine coefficients of Cl\i Cl'Zi ^3, 
<^4> bi, b 2 , b 3 , Z? 4 , and b 5 can be determined as 0.1, 6.3 x 
10" 6 , 1,995.3, 66.9, 2.34, 0.4, 2.26, 8.78, and -0.09, 
respectively. Therefore, the model for predicting the 
migration ratio of CNTs in the refrigerant-based nano- 
fluid pool boiling is expressed Eq. 20. 



£ = 1 - exp 



-66.9 



Q.lf£Y^63xl0-^-^T 
\d h J y 18/x L u b / 



fl995.3 



v 2.26" 

ds 2 PnU h \ / D \ 

184ml/ \H-L/2j 



> 3D 2 H{l-x 0 )p L _ 009 



(20) 



Figure 8a to e shows the comparison between the pre- 
dicted values of the model with the experimental data 
for different CNTs physical dimensions, refrigerant 
types, oil concentrations, heat fluxes, and liquid-level 
heights, respectively. It can be seen from Figure 8a to e 
that the migration ratio of CNTs predicted by the 
model and the experimental data have the same ten- 
dency changing with the CNTs physical dimension, 
refrigerant type, oil concentration, heat flux, or initial 
liquid-level height. The predicted values of the model 
agree with 92% of the experimental data of migration 
ratio of CNTs within a deviation of ± 20%, and the 
mean deviation is 9.96%. 

Conclusions 

Migration characteristics of CNTs from liquid phase to 
vapor phase in the pool boiling process of refrigerant- 
based nanofluid are investigated experimentally, and 
some conclusions are obtained. 

1. The migration ratio of CNTs increases with the 
increase of the outside diameter of CNTs or the length 
of CNTs. 

2. The migration ratio of carbon nanotube increases 
with the decrease of dynamic viscosity of refrigerant or 
the increase of liquid-phase density of refrigerant. Under 
the present experimental conditions, the migration ratio 
of CNTs in the R141b-based nanofluid is by maximally 
10.7% larger than that in the R113-based nanofluid, and 
is by maximally 77.4% larger than that in n-pentane- 
based nanofluid. 



Table 4 Calculation of the properties of liquid refrigerant-oil mixture 



Property 



Model for calculating property 



Author 



Specific heat 


C p ,l = (1 -x 0 )C pj + x 0 C Pi0 (A1) 


Jensen and Jackman [23] 


(J-kg" 1 -K" 1 ) 




Viscosity 


^ L _ e (x 0 ln/x 0 +(l-x 0 )ln/x r ) ( A2 ) 


Kedzierski and Kaul [24] 


(Pa-s) 




Surface tension (N-m -1 ) 


0 = o r + (a D - a r )x o 05 (A3) 


Jensen and Jackman [23] 


Density (kg-m -3 ) 
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(e) 

Figure 8 Comparison between the predicted migration ratios of the model with the experimental data, (a) for different CNTs physical 
dimensions; (b) for different refrigerant types; (c) for different oilconcentrations; (d) for different heat fluxes; (e) for different initial liquid-level 
heights. 



3. The migration ratio of CNTs decreases with the 
increase of oil concentration. Under the present experi- 
mental conditions, the migration ratio decreases by 
maximally 70.7% with the increase of oil concentration 
from 1 to 10 wt.%. 



4. The migration ratio of CNTs decreases with the 
increase of heat flux. Under the present experimental 
conditions, the migration ratio decreases by maximally 
33.9% with the increase of heat flux from 10 to 
100 kWm" 2 . 
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5. The migration ratio of CNTs increases with the 
increase of initial liquid-level height. Under the present 
experimental conditions, the migration ratio increases 
by maximally 446.9% with the increase of initial liquid- 
level height from 1.3 to 3.4 cm. 

6. A model for predicting the migration ratio of CNTs 
in the refrigerant-based nanofluid pool boiling is pro- 
posed, and the predicted values of the model can agree 
with 92% of the experimental data of migration ratio of 
CNTs within a deviation of ± 20%. 
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